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ABSTRACT 


A mathematical representation has been developed for the 
electromagnetic force field and the fluid flow field in a 
coreless induction furnace. The fluid flow field was represented 
by writing the axisvmrnetric turbulent Navier-Stokes equation, 
containing the electromagnetic body force term. The electromagnetic 
body force field was calculated by using a technique of mutual 
inductances. The k-e model was employed for evaluating the 
turbulent viscositv and the resultant differential equations were 
solved numerically. 

The theoretically predicted velocity fields were in reasonably 
good agreement with the experimental measurements reported by 
Hunt and Moore; furthermore, the agreement regarding the turbulent 
intensities was essentially quantitative. These results indicate 

t 

the k-e model does provide a good engineering representation of 
the turbulent recirculating flows occurring in induction furnaces. 

At this stage it is not clear whether the discrepancies between 
measureruents and the predictions, which were not very great in 
any case are attributable either to the model or to the measurement 
techniques employed. 
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JL INTRODUCTIO N 

In recent years there has been a growing interest in the 
develooment of an improved basic understanding of heat and fluid 
flow phenomena in coreless induction furnaces. The main motiva- 
tion for this interest is essentially twofold. One, the actual 
operation of coreless induction furnaces as moans for scrap 
melting and the refining of molten metals depends critically 
on the turbulent heat and fluid flow phenomena in this system. 

The second, more fundamental motivation is that, the turbulent 
recirculating fluid flow phenomena in this system repr-^^ent a 
very interesting class of electromagnctically driven flow problems. 

Fig. 1 shows a sketch of a typical corcless induction furnace, 
where it is seen that this consists of a refractory (or non- 
magnetic) crucible, containing a conducting metallic molt. The 
outer wall of the crucible is surrounded by water cooled induction 
coils (connected on a single pliaso or multi chase arrangement) 
through which a current is passed. The passage of this current 
through the coi) will in turn induce a current in the melt and the 
interaction of this induced current with the associated magnetic 
field will result in an electromagnetic force field, or Lorentz 
force field, which will cau.'^o in a recirculating motion of the 
melt . 

The quantitative representation of this system has two 
essential components: 

(1) The electromagnetic force field has to be calculated 
V (2) Knowing the distributed body force field the fluid 

flow field may then be obtained. 




The following general observations may be appropriate at 
this stage; 

When the magnetic Reynolds number is small the rate at which 

the electromagnetic field propagates is much faster than the 

fluid velocity thus the force field and the fluid flow field 

calculations mav be uncoupled. 

The calculation of the electromagnetic force field 

resulting from a given coil configuration is a classical problem 

in electrodynamics, which may bo readily accomplished for simple 

( 1 , 2 ) 

geometries, using analytical techniques. ' Indeed, calculations 
of this type may be readily found in the textbook literature. 

For more complex geometries, in fact when an accurate comparison 
with measurements is required, the calculations of the electro- 
magnetic force field requires numerical techniaues. 

If the flow is laminar in principle, the velocity field 

could be readily calculated by combining the known body force field 

( 3 ) 

with the laminar Navier-Stokes equations. However, this procedure 
is not readily applicable in the vast majority of practical cases, 
because the magnitude of the force field and the linear dimension 
of the system result in a turbulent flow. 

Clearly the presence of turbulence precludes the use of simple 
analytical techniques. The approaches that have been devised to 
tackle problems of this type may be divided into two major groups: 

(A) Startling with the rigorous form of the turbulent Navier- 
Stokes equations, order of magnitude approximations may be made for 
the various terms, thus approximate expressions may be deduced 
for the mean values of the various flow parameters . 
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(B) An alternative, engineering approach, which has been 
developed by Szekely and Chang, Evans and Tarapore and others^^'^ 
employed various turbulence models such as the k-W or the k-c 
model to represent the Reynolds stresses and thus obtained 
detailed maps of both the velocity fields and of the turublence 
parameters. 

Both these approaches have advantages and drawbacks. The 
technique of Hunt and Moore ' is certainly elegant and provider- 
a useful insight into relationships between the key system 
parameters. However, because of the approximations involved the 
model cannot be used to predict the detailed velocity field, 
the maps of the turbulence and thus cannot but provide very pre- 
liminary order of magniture estinutes regarding the key heat, and mass 
transfer and dispersion phenomena that are of primary interest in the 
solution of practical engineering problems. The numerical approaches 
that have been developed for tackling problems of this type are 
certainly attractive because they are capable of addressing the 
very questions of practical interest, such as local shear rates, 

dispersion rates and ultimately the transport controlled reaction 

.. ( 8 , 12 - 15 ) 

kinetics. ' • 

The main drawbacks of these engineering calculations include 
the substantial computational labor required and the lack of a 
fundamental basis for the turbulence models employed. 

In view of the great potential usefulness of these techniques, 
from a pragmatic standpoint it would be highly desirable to develop an 
a posterioiri justification for their use, through a critical 
comparison of the theoretical predictions with experimental 


measurements . 
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Up to the present only a limited range of such measurements is 
available, which encompases data obtained both from laboratory 
scale apparatus and industrial scale units. Many of these 
measurements involved the determination of surface velocities 
and tracer dispersion rates; while the agreement between oredictions 
and the measured data was quite good, Auguring well for many j 
engineering applications a really rigorous test of these models 
has not been possible because of the lack of sufficientlv accurate 
measurements. 

In a recently published paper, which was kindly provided 

for us prior to formal publication. Hunt and Moore have reported 

on an interesting set of measurements in an inductively stirred 
(6 7) 

mercury system ^ ' which should provide a useful basis for a more 
rigorous test of the model. The purpose of the paper is to report 
on such an assessment. 

2. EXPERIMENTAL 

The experimental measurements, which will be used for the 
purpose of comparison have been reported by Hunt and Moore and 
therefore only a very brief recapitulation will be given here. 

A schematic sketch of the experimental arrangement is 
shown in Fig. 2. In essence the apparatus consisted of a water 

cooled stainless steel vessel annroximatel'' 0.3m in diameter and 0.4 m 
high, containing mercury. Agitation was provided by an induction 
coil containing 11 uniformly spaced turns. The coil current was 1900 A, 
RMS, having 50 cycles. 

The actual measurements taken included the determination of the 
electromagnetic force field, using search coils and the measurement 
of both the time smoothed and the fluctuating velocities, using a 
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mechanical probe. In essence this mechanical probe consisted 
of a perforated spherical shell, made of tantalum and the actual 
measurement involved the determination of the drag exerted on this 
sphere. Measurements close to tank wall were made with "wall 
probe" described elsewhere. . 

The actual results of the experimental measurements will be 
discussed in a subsequent section. 

3. THE ANALYSIS 

The salient features of the analysis will be given in the 
following: 

As mentioned earlier the analytical oroblem is to calculate 
the electromagnetic force field and the turbulent fluid flow field 
in i- cylindrical system, agitated by a symmetrically placed induc- 
tion coil. 

3.1 Calculation of the Electromagnetic Force Field 

The electromagnetic force field is given by: 

F = J X B (1) 

where 

J is the induced current density and 
B is the magnetic flux density 

Thus the problem is the calculated J and B for a given geometry 
and coil current. 

For idealized systems, e.g. infinitely long cylinders, and a 
travelling wave, this task may be accomplished analvtically , 
through the solution of Maxwell’s equations However, for the 
present system, particularly when a quantitative comoarison is 
desired, an alternative technique, employing the concept of mutual 
inductances is preferable. 
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For a cylindrical crucible which is characterized by axial 
syminetry, the lines of equal current density are circles in 
planes perpendicular to the axis of the coil» i.e. the 6 direc- 
tion. From Amper's law the magnetic field in an elementary 
circuit in the melt may be written in terms of the vector 
potential A(B=VxA) as; 






dH coil 

d — ^ + r 
^ k=l 


I(k) rf 



(2) 


where ^gis the vector potential, d£ is the line element of a 
circuit of constant current density and S is the cross sectional 
area of the circuit. I (k) is the coil current. 

It should be noted that first term on the rhs of Eq, (2) 
describes the induced potential while the second term described 
the applied potential. 

From Faraday's law, the induced current in this circuit 
can be written for a time harmonic field as; 


where 


J. 

X 
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is the mutual ino.uctanre 


(3) 

(4) 


The calculation then ••>rocccds by dividing the r.clt into elementary 
circuits (14x14 in the present case), each being represented by Eq. (3) 
to be solved simultaneously, to obtain the induced current. Once the 
induced current field is known the vector potential may be obtained 
from Eq. (2) , when then facilitates the computation of the electro- 
roa'‘-">,c'c,ic hod*> force field. 
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3.2 Calculation of the Fluid Flow Field 
For cylindrical symipetry, the equation of continuity and 
motion take the following form: 

equation of continuity: 


-^iprUd+T^U-)=0 ‘ 

r cr cz 

Momentum balance in the z-direction 




Re(JgXBp/2 


Momentum balance in the r-direction 


(5) 


( 6 ) 




(7) 




Re(JgXBJ) 


Here the quanity is the effective viscosity defined as: 


^eff 


y + y^ 




where 

= (9) 

is the turbulent viscosity. 

Here k and z are the turbulent kinetic energy and the turbulent 
kinetic energy dissioation respectively. 

Separate transport equations have to be written down for these, 
which take the follov;ing form: 


ORIGINAL PAGL tS 
OF POOR QUALITY 
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.( 10 ) 


where S(^ is the net rate of generation of turbulent properties 
per units volume. The source terms Sj^and of transport equation 
of k and c are given by 

Sj^ = G - D .(11) 

= Cl g/k G-C 2 P e^Vk (12) 

where 


3U 

G=„^12((_ 


D = p e 


3U 


+ ij^)^ + (-4-)^{+f 


U 




3U 


~+ -^)^] 
3r dz ‘ 


.( 1 . 3 ) 

(14) 


The value used for the constants, viz C^(=0.09), C^(=1.45) 

and C2(=1.92) were taken from the work of Spaulding and Launder, 

(13,17) and v.'ere not adjusted in any wav in course of the computation. 


The boundary conditions, which are required to complete 
the statement of the problem will have -O specify zero velocity 
at the solid surfaces, zero shear at the free surface and the 
existence of symmetry about the centerline. Regarding the 
quantities k and c these were aet equal zero at solid surfaces 
and their gradient was stipulated to be zero at the free, surfaces 
It is noted, furthermore that wall functions were used to 
define the shear stress at the solid surface. These wall 

functions were based on universial velocity distribution for 
turbulent boundary laye.rs without pressure gradient and body forcer.. 
This is an area where careful experimental measurements could 
provide a useful refinement. 

The governing equations were put in a finite difference 
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form, using a 16x16 grid and the resultant set of simultaneous 
non-linear algebraic equations was solved using an iterative 

technique. Typically the computer time requirements were about 
120 seconds on MIT's IBM 370 digital computer. 

4. COMPARISON OF THE COMPUTED RESULTS WITH MEASUREMENTS 

Since the principal concern here is a critical^ comparison 
of the experimentally measured and the theoretically predicted 
velocity fields and turbulence characteristics the treatment 
presented will be confined to this aspect of the problem. In any 
case tiie calculation of an electromagnetic force field for a given 
coil configuration may be rerdily undertaken as an essentially 
routine matter. 

Figs. 3 and 4 show the experimentally measured and the 
theoretically predicted maps of the velocity vector. Inspection 
of these two figures shows quantitative agreement regarding the 
nature of the flow and the position of the recirculating loops. 

While this form of representation is not ideal for making a fully 
quantitative comparison between the numerical values of the velocity 
vector, these magnitudes seem comparable. 

Fig. 5 shows the experimentally measured profile of the 
axial velocity, determined at a vertical position, corresponding 
to the eye of the upper vortex. Also shown, with the broken 
line is the computed velocity profile. It is seen that these 
two profiles are quite similar, but that the two curves do not 
coincide, vrhether this discrepancy is attributable either to 
possible shortcomings of the model or to experimental errors will 
be discussed subsequently 
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Fig. 6 shows a plot of the local value of the turbulence 
intensity, as a function of the radial position. Here again 
the measurements are given with the circles, while the 
theoretical predictions are designated by the full line. It is. of 
interest to note that here the agreement is very good, rather 
better than that found for the velocity profiles. 

5. DISCUSSION 

In the paper a comparison is made between experimentally 
measured fluid velocity profiles and profiles of the turbulence 
intensity as reported by Hunt and Moore for an inductively stirred 
mercury pool, and theoretical predictions based on the concept 
of mutual inductances and the two equation k-e model. 

This comparison was thought to be instructive, because 
these measurements, kindly supplied by Messrs H nt and Moore 
represent the most detailed data on inductively stirred systems 
available up to the present. 

The principal findings of the work may be summarized as 
follows : 

(1) The technique of mutual inductances was able to 
predict the electromagnetic force field quite readily, thus this 
facet of the approach was unlikely to introduce a serious error. 

(2) The theoretical predictions regarding the general 
nature of the flow were found to be in very good agreement 
with the measurements. 
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(3) The profiles of the time smoothed velocity near the 
wall were reasonably well predicted by the model although there 

were certain discrepancies, oatricularly at large distance from the 
wall . 

(4) The agreement between measurements and oredictions, 
regarding the profiles of the turbulence intensities was very good. 

Before considering these points in detail some general 
observations may be appropriate. 

As far as the mathematical modelling of turbulent 
recirculating flows is concerned it is fair to say that important 
reservations must be expressed regarding the fundamental basis 
of the k-£ model and that of the wall functions. 

It should be stressed, however, that at present there 
appears to be no generally applicable alternative for doing 
engineering type calculations for turbulent recirculating flows 
where details are required of the turbulence energy distribution 
and of tho local heat or mass transfer rates. 

Thus in a pragmatic sense the application of this mpdel 
has to be justified through a direct comparison with measure- 
ments. On the basis of information available in the literature 
the use of this model has been quite successful in a number of 
instances. 

When applying turbulent recirculating flow models ‘to 
representing electromagnetically driven flows an - Mitional 
complication may arise because of the possible damping or 
accentuation cf anistropv^as aused by the interaction between 
the electromagnetic force field and the turbulence field. 
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On the basis of prior experience one would expect 
turbulence models of the type employed here to predict the 
overall flow field rather well, because neither the particular 
features of the turbulence model employed nor the damping effect 
of the electromagnetic field on the turbulence are likely to be 
very important in the bulk. More specifically ir. the bulk of the 
fluid the convective transport of momentum is likely to dominate. 

It follows that a much more critical test of the model would be 
provided by the assessment of the predictions in the near wall 
region. 

In considering the experimental technique employed, it 
has to be recognized that the characterization of electromagnetically 
driven turbulent recirculating flows is notiorously difficult, 
because of the problems inherent in the use of traditional measuring 
equipment. The technique employed by Moore and Hunt is thought 
to be ingeneous, but perhaps not subjected to a lengthy enough 
testing period to eliminate all possible experimental errors. The 
relatively large size of the probe or even the alternative prob- 
ing devices employed near the wall compared to a standard hot 
film device) makes the reliability of the measurements in the very 
near wall regions somewhat problematic, where accurate data would 
be most desirable. 

It follows that it is not really quite clear whether the 
discrepancies between measurements and predictions are unequi- 
vocally attributable to the shortcomings of the model or to the 
measurement techniques. 





! 
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Let us now comment briefly on the specific findings of this 
study, concerned with the details of the velocity fields and 
of the spatial distx'ibut ion of the turbulence intensity. 

While there was good agreement between the predicted and 
the experimentally measured velocity profiles a certain discrepancy 
was evident, both in the immediate vicinity of the wall and in the 
bulk of the fluid. 

It should be remarked that the experimentally obtained 
velocity profile do not satisfy the equation of continuity i.e. 

R 

^ U^2Tirer = 0 (15) 

thus some questions may be raised regarding the accuracy of 
the absolute values of the measured velocities. At this stage 

one may comment that tlie calibration curves for the price were 
not exactly linear, which in conjunction with the straight line 
relationship actually employed provided a significant error band 
especially at low velocities. It should be stressed here that this 
should not be taken as a criticism of the experimental technique 
because the problems of velocity measurements in liquid metal systems 
are generally appreciated by workers in this field. 

It is of interest to note that perhaps unexpectedly very good 
agreement between the theoretically predicted and the experimentally 
measured turbulence intensities. The k-t model explicitly assumes 
local anisotrophy and serious questions may be raised, wl^other 
this condition has been met in an inductively stirred system. One 
may speculate that some self-cancelling errors may have come into 
play, alternatively perhaps the experimental technique was not 
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sensitivc enough to pick up some of the higher harmonics. It should 
be remarked, furthermore, that the errors introduced due to the non- 
linearity of the calibration curve will affect the absolute values 
of the velocity far more than the velocity ratios. Since the 
turbulence intensity curves, shown in Fig. 6 represent ratios, these 
would not be less affected by calibration errors. 

In conclusion one must state that notwithstanding its short- 
comings as far as lacking a fundamental basis, the k-c model 
appears to provide a reasonable prediction for both the gross 
features and the detailed velocity fields and turbulence energy 
distribution in an inductiv'cly stirred system. 

The principal discrepanev that o.ie should expect would lie 
in the near wall regions. Should one require more precise 
information regarding the behavior of these regions such as the 
knowledge of local heat or mass transfer rates then the best 
alternative would be to measure these quantities directly. Work 
of this type is in progress in the author's laboratory at present. 
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Syml )0.1j> 

A ’ vector potential f w/b/m) 

2 

B magnetic flux density (W b/m ) 

electromagnetic stirring force (N/m^) 

I coil current (A) 

2 

J ■ induced eddy current (A/in ) 

2 2 

K turbulent kinetic energy (in s ) 

2 

p pressure (N/ni ) 

K radial coordinate • . 

. 2 
S cross section area of elementary circuit (m ) 

z axial coordinate 

magnetic pormcabilitv (H/m) 

II effective viscosity (kg/m.s) 
off 

p density of the moll (kq/m'^) 

2 3 

f turbulent energy dissipation (m s ) 

o electric conductivity (Ohm.m) ^ 
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Figure Capti ons 

1. Sketch of an induction furnace 

2. Salient features of the apparatus 

3. Experimentally measured velocity field 1900 Am;j. coil current 

4. Computed velocity field for 1900 Amp coil current 

5. A comparison of the experimentally measured (full lino) 
and theoretically r^redicted broken line) radial variat.ion - 
of tlie axial velocity component at a height of 22.5 cm. 

6. A comparison of the experimentally measured (circles) and 
the theoretically predicted (full line) radial profile 

of the turbulence intensity. 









ORIGINAL PAGE IS 
OF POOR QUALITY 









Fiq. 2 









30 


25 

20 


15 



5 


ORIGINAL PAGE IS 
OF POOR OUAUTY 


t \ f t 7 

i t f r / / 

] 1 1 1 /////• ^ 

Hi., (///// 

\ \ \ \ \ I ; 



n n 


1 






\ I ^ ^ 

1 1 \ \ \\ \ 

^ 1 M \ 


\ 

\ 

^ \ 

^ I 

1 

y 

^ / 

' / 


\' 


\' 


\ \ N \ \ \ 


// 






\ t •e' ^ 



\ \ J / / /■"/// /^ 

i i / / / I //// /^ 

M 1 / / \ ui / I / 


n / i / / \ 1 1 1 1 ' 

I t I I I 1 * • » I « 


1 


I 


/ 


f 

i 


i I 


i 


I 


i 


1 1 1 

, I ^ ^ 

I \ \ ^ ^ 

\ \ \ \ \ ^ '■ 

M \ \ \ \ \ '" '* 

1 \ 

» I VI \ '■ 

t I V ‘ 

V \. \ 











I 


— measured 

— predicted 





ORIGINAL PAGE IS 
OF POOR QUALITY 


— predicted 
0 measured 



L 


0 


3 


J L_ 

6 9 

r cm 


Fig. 6 


/ 




5 


r,'- f. O- 

I.' ‘ ^ 




; I • 

‘ : . t • 



